arXiv: 1505.05092v2 [hep-th] 15 Sep 2015 


SU-ITP-15/12 


Prepared for submission to JHEP 


Dissecting holographic conductivities 


Richard A. Davison^ and Blaise Gouteraux^’^ 

^Lorentz Institute for Theoretical Physics, Niels Bohrweg 2, Leiden NL-2333 CA, The Nether¬ 
lands 

‘^Stanford Institute for Theoretical Physics, Department of Physics, Stanford University, 
Stanford, CA 94305, USA 

^APC, Universite Paris 1, CNRS, CEA, Ohservatoire de Paris, Sorbonne Paris Cite, F- 
75205, Paris Cedex 13, France 

E-mail: davison@lorentz.leidenuniv.nl, gouterau@stanford.edu 

Abstract: The DC thermoelectric conductivities of holographic systems in which 
translational symmetry is broken can be efficiently computed in terms of the near¬ 
horizon data of the dual black hole. By calculating the frequency dependent con¬ 
ductivities to the first subleading order in the momentum relaxation rate, we give a 
physical explanation for these conductivities in the simplest such example, in the limit 
of slow momentum relaxation. Specifically, we decompose each conductivity into the 
sum of a coherent contribution due to momentum relaxation and an incoherent contri¬ 
bution, due to intrinsic current relaxation. This decomposition is different from those 
previously proposed, and is consistent with the known hydrodynamic properties in the 
translationally invariant limit. This is the hrst step towards constructing a consistent 
theory of charged hydrodynamics with slow momentum relaxation. 
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1 Introduction 

Understanding the transport properties of strongly interacting many-body systems with 
no quasiparticles is a topic of much interest for both experimental and theoretical 
reasons. One class of theoretical examples are the strongly interacting quantum held 
theories which are holographically dual to classical theories of gravity. Holographic 
duality can be exploited to calculate the transport properties of these examples in a 
relatively simple way, with the goal of determining non-holographic effective theories 
which control these properties. Recent examples of this approach include [1-9] . In this 
paper, we study holographic systems with weakly broken translational invariance as a 
first step in formulating a general hydrodynamic theory of strongly interacting systems 
with slow momentum relaxation. 

The transport properties of primary interest are the electrical (cr), thermoelectric 
{a) and thermal (k) conductivities that control the linear response of the electric current 
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J and the heat current Q to small electric helds E and temperature gradients VT 


( J\ _ ( (T aT\ ( E 
Vq; \aT rt) V-vt/t 


( 1 . 1 ) 


The primary consideration in determining the qualitative form of these conductivities in 
holographic systems is whether the total momentum P of the system is approximately 
conserved or not. In this paper, we will primarily address situations in which this is 
the only long-lived quantity. In these cases, a perturbative expansion in the (small) 
momentum relaxation rate T can be performed within the memory matrix formalism 
[1, 2, 8, 10, 11]. To leading order in this expansion, the conductivities are all Drude-like, 
with DC values determined by T and by the static susceptibilities xjp and xqp of fho 
translationally invariant state where momentum is exactly conserved: 


a{u) = 


Xjp 


xpp r 


lU} 


/ N XjpXqp 
a (cu) = ——^ — 

Txpp r 


lU} 


K (cj) = 


X-%p 


Txpp r 


lUl 


( 1 . 2 ) 


Physically, any current A which overlaps with the momentum {xap 7^ 0) cannot decay 
at a rate larger than T at late times. The slow relaxation of momentum acts as a 
bottleneck that forces the current into a coherent late time response, even if its intrinsic 
relaxation timescale is fast. The low energy optical conductivity is dominated by a 
single pole that is parametrically close to the origin. In the opposite situation, when 
a current does not overlap with the momentum {xap = 0), it will dissipate at its 
intrinsic rate. This is an example of incoherent transport [6] and is the case for the 
electric conductivity in charge conjugation symmetric states, for example. 

These results form a basic, non-holographic effective theory that describes the 
transport of charge and energy in holographic systems in which momentum is approx¬ 
imately conserved. For these systems, one can attempt to enhance this basic effective 
theory by combining it with our knowledge of the system’s properties in the transla¬ 
tionally invariant limit, in which its late time behaviour is described by the laws of 
relativistic (conformal) hydrodynamics. The simplest way to incorporate the above 
results from the memory matrix formalism is to modify the momentum conservation 
equation in hydrodynamics, such that P decays at a constant rate T. This yields the 
conductivities [1] 
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where e, p, n, /i and s are the energy density, pressure, charge density, chemical poten¬ 
tial, and entropy density of the state respectively. Each conductivity has a coherent 
contribution at leading order in E, as well as a subleading incoherent contribution pro¬ 
portional to the intrinsic conductivity of the hydrodynamic state aq.^ The former is in 
perfect agreement with the memory matrix results (1.2), while the latter is a correction 
due to long-lived diffusive modes, whose form is specihed precisely by the relativistic 
hydrodynamic theory in terms of a single transport coefficient aq. The memory matrix 
results (1.2) can be extended to incorporate the effects of diffusion in more general 
setups [8] ~ these are independent of the momentum relaxation rate and enter at the 
hrst subleading order in a small T expansion. 

Recent advances in the holographic description of strongly interacting systems with 
momentum relaxation, in particular, the discovery of analytically tractable toy models 
of such systems [12-14], and the development of efficient calculational tools to determine 
their DC conductivities [15-19], have made it easy to test this modihed version of 
hydrodynamics. At leading order in F [7, 20-22] , the holographic results are consistent 
with those of the memory matrix (1.2) and therefore with the leading order modihed 
hydrodynamic results (1.3). However, the holographic results are inconsistent with the 
modihed hydrodynamic results at sub leading order. To be explicit, we will consider 
the gravitational action [14] 


S = J d^x^g (n + 6- ^ ^ dyjd,cl>j^ , 


which has the analytic black brane solution [23] 
ds^ = -r^f{r)dt‘^ P {dx‘^ + dy^) , 

j- 3 2rQ 4rQ / 4r‘^ 


01 = mx, 02 = friy, 
At{r) = 


(1.4) 


(1.5) 


The massless scalar helds 0/ break translational symmetry and so the DC conductivities 
are hnite [14, 18] 
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where the explicit expressions for the energy, charge and entropy density of this state 
are given in (2.7) and (2.8). The parameter m controls the strength of translational 

^Note that ctq is an intrinsic property of the finite density state obtained by perturbing the neutral 
UV CFT by a chemical potential. 
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symmetry breaking and therefore the rate of momentnm relaxation in the dnal field 
theory state. 

In [5, 15, 17-19], it was suggested that anc could be interpreted as being composed 
of two physically distinct pieces: a coherent contribution /m? due to momentum re¬ 
laxation, and an incoherent contribution 1 (see [24, 25] for further related work on the 
frequency dependence of the thermoelectric conductivities). However, this is incon¬ 
sistent with the known value of the incoherent contribution ctq in the translationally 
invariant (m = 0) limit [26] 


(1.7) 

m=0 

Furthermore, it is clear that a decomposition of this kind is inconsistent with the other 
hydrodynamic DC conductivities, as can be seen by comparing (1.6) with the DC limit 
of (1.3). 

In this paper we resolve these problems, and for the first time provide a clear 
description of the physical processes underlying the simple DC conductivities (1.6), 
by analytically calculating the low frequency conductivities for the holographic theory 
(1.5) at small values of m where there is approximate momentum conservation. We 
identify two physically distinct contributions to each conductivity - a coherent contri¬ 
bution controlled by the slow relaxation of momentum, and an incoherent contribution 
due to the intrinsic conductivity uq. The value of ag we obtain is consistent with the 
known value in the translationally invariant limit (1.7) [26]. Technically, we achieve 
this decomposition by changing basis from the currents J and Q to more theoretically 
convenient currents J± which are orthogonal: the two-point retarded Green’s function 
of J+ with J_ vanishes. To the first subleading order at small T the conductivity of 
J_ is entirely coherent, and that of J_|_ is entirely incoherent. In the strict T = 0 
limit, these reduce to the currents P and 3eoJ/2 — tt-qP respectively, which decouple 
and capture the entirely coherent and entirely incoherent responses respectively in the 
translationally invariant, hydrodynamic system (subscript Os here denote the thermo¬ 
dynamic quantities of the m = 0 state). 

To the hrst subleading order at small u and T, with ca/T hxed, we find that the 
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frequency dependent conductivities take the form 
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where the momentum relaxation rate is 

r = —r (1 + Xm? + 0{rri^)) , 

47r(e + p)^ ' 


(1.9) 


the thermodynamic quantities are those of the m 7 ^ 0 state, and A is given in equation 
(3.11). For comparison with the hydrodynamic results (1.3), these expressions may be 
written as 
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( 1 . 10 ) 

( 1 . 11 ) 

( 1 . 12 ) 


These results isolate the reason for the inconsistency between the modihed version 
of hydrodynamics and the holographic system: the modihed version of hydrodynamics 
does not adequately describe the coherent component of the system’s response. Al¬ 
though it reproduces the correct coherent contribution at leading order in F, it does 
not adequately account for the hrst subleading corrections to this. These corrections 
are important as they enter at the same order as the incoherent contribution, and em¬ 
phasize the need for a more systematic derivation of how hydrodynamics is modihed 
by the weak breaking of translational symmetry. Our calculation also highlights the 
important message that it is in general not possible to separate the coherent and inco¬ 
herent contributions to the conductivities from their DC expressions (1.6) alone. We 
note that the obvious decomposition of the DC conductivities (1.6) still has physical 
meaning in terms of the DC conductivities at zero electric or heat current [11, 18]. 

Finally, although the main focus of our paper is the limit of slow momentum 
relaxation, we can also easily access the regime of fast momentum relaxation in the 
holographic theory (1.4). This is a regime in which neither the hydrodynamic (1.3) 
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nor memory matrix results (1.2) are applicable. For any value of m, it is possible to 
diagonalise the response of the currents by changing to an appropriate basis J±. In the 
limit m —)■ cxD, the decoupled currents are precisely J and Q, the electrical and heat 
currents. It would be very interesting to determine whether a low energy decoupling 
of this type is present more generally in systems with fast momentum relaxation (in 
particular, those with a potential for the scalar helds, which are more reliable from the 
point of view of string theory). 

In Section 2, we identify the diagonal J± basis of currents in the field theory by 
decoupling the bulk field equations, and examine how these decoupled currents relate 
to J and Q in various limits of interest. In Section 3, we determine the frequency 
dependence of the conductivities of the currents J± in the limit of slow momentum 
relaxation, showing that one is coherent and that one is incoherent, and explain what 
this means for the conductivities of J and Q. We conclude in Section 4 with an outlook 
for future work. The appendices contain some technical details of our holographic 
Green’s function calculations. 

2 Diagonalisation of the conductivities 

To determine the frequency dependent thermoelectric conductivities in the strongly 
interacting field theory state dual to (1.5), we will use the Kubo formulae [1] which 
relate these conductivities to the retarded two-point functions of the currents J 
and Q = Je — ^ where Je is the energy current: 

cr (ca) = — [Gjj (cj, fc = 0) — G^j (ca = 0, /c —)■ 0)] , 

OJ 

“ (i^) = ^ [G'qj k = Qi) — Gqj {uj = 0,k ^ 0)] , (2.1) 

k (i^) = ^ [^QQ k = 0) — Gqq (o; = 0, fc —)■ 0)] . 

To evaluate the Green’s functions on the right hand side, it is convenient to first change 

the basis of currents and not work directly with J and Q, as we will shortly describe. 

Following that, we will use the standard tools of the AdS/GFT correspondence to 
compute the Green’s functions. 
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2.1 Decoupling of the gravitational equations of motion 

To determine the two-point functions of J and Q, we consider the following consistent 
set of linear perturbations around the black brane solution (1.5) 

'59?(*.r)= [ , Sgt{t,r) ^ , 

rZ 

ay{r,uj)e 6(()2{r,t) = j —X2{r,uj)e 

where indices are raised with the background metric. These obey the following lin¬ 
earised equations of motion (where primes denote derivatives with respect to r) 
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(2.3) 


2 / 
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which are comprised of two linearly independent dynamical equations, and one con¬ 
straint equation. 

We can decouple the two dynamical equations by changing variables to 


where 


Here 


^ [r‘^fZ'±\ + ^ + r^ff + m-f±A[rf) = 0, 


'ip± = — {hi + iuhl) + rA!^ay\ + '^±ay . 
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with e and n, the energy and charge densities of the state, given by 

n = fiVQ. 


j - 2rn 1 - — + — 

" ' 2ri ^ 4r?, ' ’ 


(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 


0 ^'0, 

The other thermodynamic properties of the equilibrium state are [14] 

TTl^ ■ ■ 1 

dvrT = 3ro-, s = Anri , p = (T“) -|- r^mZ = -e + r^mZ. (2.8) 

2ro 4ro 2 
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This change of variables corresponds, in the field theory, to a change of operator 
basis^ from to (J+, J_). The deconpling of these variables in the bnlk corre¬ 

sponds to a diagonalisation of the matrix of two-point fnnctions of the dnal operators 
i.e. it corresponds to diagonalising the matrix of condnctivities. By a carefnl analysis 
of the on-shell action (see appendix A), we hnd that the two-point retarded Green’s 
fnnctions are related to the bonndary behavionr of the deconpled helds 'ijj± via 

{JJ){uj) = -^— ( 7 + 0 + M - 7-0- M), 

7 + - 7 - 

7Tl 

{JJe){u}) = {JeJ){u}) = -(0+ (w) - 0_ (w)) + n, (2.9) 

7 + - 7- ^ ^ 

771 ^ S 

{JeJe){u:) = -(7+0- { 00 ) - 7-0+ M) + x6, 

7 + - 7 - 2 

where the angled brackets denote the specihc combinations of retarded two-point Green’s 
fnnctions which enter in the Knbo formnlae for conductivities (2.1) 

{00){u) = Ggc,(a;, k = 0) - = 0,k^ 0), (2.10) 

and where 

0 + M = , (2.11) 

W± 

r^oo 

are determined by solving the decoupled equations of motion with ingoing boundary 
conditions at the horizon, and contain all of the frequency dependence of the correlators. 

It is simple to diagonalise the matrix of correlators by dehning the currents (with 
overall normalisation constants a±) 

J± = a±{JE+ ^±^17), (2.12) 

so that {J±J±) depends only on 0±, and the cross correlator {J±Jzf) = 0. Physically, 
this diagonalisation of the matrix of conductivities means we can divide the response of 
the currents of our system, at any energy scale, into two completely independent sectors, 
each with its own spectrum of excitations. This situation is familiar, for example, 
in zero density, translationally invariant systems, where J and Je decouple due to 
charge conjugation symmetry. In our case, there does not appear to be any symmetry 
protecting this exact decoupling at all energy scales, and we do not expect it to be 
true in general for holographic systems. The more pertinent question is whether the 
conductivity matrix can be diagonalised at low energies u in more general holographic 
states. This does not necessarily require an exact decoupling of the bulk perturbations. 


^The energy current in our system is Je = = 0), which is the momentum when m = 0. 







and would be an indicator of the existence of a simple, low energy effective description 
of transport in these states. An example of this, when m = 0, is described below. 

The perturbation equations for this holographic model can also be completely de¬ 
coupled at non-zero wavevectors k. Again, we expect that this feature is specihc to 
this very simple example, and will not be true in general. 

Inverting these relationships, we can express the responses of the correlators we are 
truly interested in - those of the electrical and heat currents - as linear combinations 
of those of the decoupled currents J± as follows: 






m i 
2 




{QJ) = {JQ) = -7- + - {JM + 7+ + - {J-J-) 


m 
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{QQ}^(i- + -) {J-+J+) + 7++-) {J-J-), 

V m/ V m/ 

where we have introduced the rescaled correlators 




=pm 7± 


al (7+ - 7_)" 7+ - 7- 


0± (w) + 


(7+ - l-Y 


(2.13) 


3 ^ 

+ 2m7±n ) , (2.14) 


for convenience. From these, one simply needs to divide by the appropriate factor in 
the Kubo formulae (2.1) to extract the relevant conductivity. The decoupled currents 
J'± are sourced by ± {E/m + (yzp -|- ^/m)'VT/T), and transport the conserved charge 
densities (T** -|- 7 ±mJ*) / ( 7 + — 7 _). 


2.2 The decoupled currents in various limits 

Remarkably, we have managed to decouple the response of the currents J± at all fre¬ 
quencies, and for all values of the parameters m, T and fi. For certain values of the 
parameters, the decoupled currents J± take particularly simple forms, which allows us 
to ascribe a clear physical meaning to the decoupling. 

The simplest limit is already very familiar: when /i —)■ 0 at fixed T and m, after 
an appropriate choice of normalisations the decoupled currents are J+ —)■ J and J_ —)• 
Je = Q- This is simply the well-known result that at zero chemical potential, the heat 
and charge currents of a system decouple due to charge conjugation symmetry. The 
heat and charge conductivities may be qualitatively different from each other in this 
limit. The charge response will be incoherent as J does not overlap with any almost 
conserved operators. As Q overlaps with P, the heat response will be coherent when 
P dissipates slowly (at small m), and incoherent otherwise. 

There is another limit which is in fact rather similar to this: when m —)■ 00 at hxed 
T and p, after an appropriate choice of normalisations, the decoupled currents are just 
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J+ —)■ J and ^ Je — = Q- This is a rather surprising result: in the limit of 

very strong translational symmetry breaking, the charge and heat currents decouple! 
Heuristically, it is as if there is an emergent form of charge conjugation symmetry in 
this limit. One way of understanding this is that when m ^ /i, T, the contributions 
of the uncharged scalar operators dual to the helds 0/ dominate the thermodynamic 
properties of the system such that it looks like a neutral state. In particular, the 
dimensionless ratio of charge density to entropy density, a thermodynamic measure 
of the ratio of charged to neutral degrees of freedom, approaches zero in this limit; 
n/ s ~ /i/m —)■ 0. However, this is qualitatively different from the /i = 0 limit in that it 
is specihcally the heat current Q which decouples from J, while other neutral currents 
like Je still couple to J. It is clearly worth investigating to see if this a common feature 
of low energy transport in states of this type, or just a peculiarity of this holographic 
system. Finally, note that in contrast to the previous /i —)■ 0 limit, in this limit both 
the charge and heat conductivities will be incoherent, as momentum dissipates quickly 
in the limit m —)■ cxd. 

Finally, there is the limit of slow momentum relaxation, in which we are mainly 
interested in the remainder of this paper. In the limit m —)■ 0 with T and /i hxed, the 
decoupled currents asymptote to 

Q 

t/-|_ — ^ Je — t;— J + 0(jn^), —y Je ~\~ O(m^), (2.15) 

2no 

after appropriate choices of normalisation, where the subscript Os denote the thermo¬ 
dynamic quantity in the m = 0 state. This decoupling also has a clear physical origin, 
which is independent of holographic duality. In the strict m = 0 limit, our state obeys 
the laws of conformal, relativistic hydrodynamics at low energies. In such a hydrody¬ 
namic state, the currents given in (2.15) decouple at low energies. Je = P controls 
the coherent component of the system’s response while the other current controls the 
incoherent response, since it decouples from the conserved total momentum P. See [27] 
for more details. 


2.3 DC contributions of each sector 


The DC conductivities correspond to the a; —)■ 0 limits of the subtracted correlators, 
and are given analytically in (1.6). We will conhrm these results in the next section. 
From these, we can extract the DC limits of the diagonal correlators 


Im 


lim -{J±J±) (cj) 

uj ->-0 00 




2n 


3e./l + 


9e2 


(m^ + /i^) i 



+ AirnT 


(2.16) 
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Using the decomposition (2.13), we can then easily isolate how much each of the decou¬ 
pled sectors contributes to each DC thermoelectric conductivity. Although the full DC 
conductivities are very simple, each individual contribution is given by a very compli¬ 
cated expression (which can be found by combining (2.13) with (2.16)). The separation 
of the full conductivities into two decoupled sectors is highly non-trivial and cannot be 
guessed just from the form of the DC conductivities. 

It is instructive to examine these contributions in the various limits of the previous 
subsection. We use the notation that the superscript ± indicates the contribution of 
the J± sector to each conductivity. In the p —)• 0 limit (at hxed T,m), 


^DC ^DC + 

a DC ), 


^DC 

“dC 

(2.17) 

^DC 

Rdc ^DC + 



This limit is well-known and the the nature of the decomposition is clear: this is 
the charge conjugation symmetric limit in which the decoupled currents J+ and J_ 
are the charge and heat currents respectively. Therefore, the electrical and thermal 
conductivities are controlled completely by the -|- and — sectors respectively, and the 
off-diagonal conductivity vanishes, to leading order at small p. 

From this point of view, a qualitatively similar limit is the limit m —)■ cx) (at hxed 
T, /i), where 


^DC + 

^DC 0 {m-^) , 


(^DC 


(2.18) 

^DC 0{m-^), 

^DC ^DC + 0{m~^). 



Again, this is easy to understand: the decoupled currents in this limit are again the 
charge (J+) and heat (J_) currents, and so the electrical and thermal conductivities 
are hnite^ and determined, at leading order, by the -|- and — sectors respectively, while 
the off-diagonal conductivities vanish at leading order in this limit. 

Finally, let us turn to the limit we will address in the remainder of the paper: 
m —)■ 0 (at hxed T,/i). This is the limit of slow momentum relaxation. In this limit, 
the contribution of each of the sectors to the DC conductivities take a very suggestive 

^It is crucial here that the horizon radius is replaced by its expression in terms of physical parameters 
T, /i and m before taking the m —>■ oo limit. When this is done, R, does not vanish at large m, contrarily 
to what the formula in (1.6) might appear to indicate. 


11 



form 



^DC ^DC - 0-Q + O(m^), 


DC 

«DC + O(m^), 

(2.19) 


Kdc I^DC - + 0{rn^), 



where aq is given in equation (1.7). Each DC conductivity clearly decomposes into two 
independent pieces, one of which is exactly equal (to this order in m?) to the incoherent 
contribution, due to intrinsic current relaxation, present in the effective hydrodynamic 
theory (1.3). Although the full holographic DC conductivities are not consistent with 
this effective hydrodynamic theory, the decomposition above suggests that at least this 
part of the hydrodynamic theory is accurate. To conhrm this, and give a more physical 
interpretation to the decoupling in this limit, we will now calculate the frequency 
dependence of each of the two independent contributions to the conductivities. 


3 Frequency dependent conductivities at small m 

To determine the frequency dependent conductivities, we need to solve the perturba¬ 
tion equations (2.4) at non-zero u, with ingoing boundary conditions at the black brane 
horizon. We can only End analytic solutions to these equations by working perturba- 
tively in a small frequency expansion. The resulting conductivities we hnd, extracted 
via (2.13), are fractions, with both the numerator and denominator given by power 
series in cu. This procedure was used in [28-30] to determine the two-point Green’s 
functions of translationally invariant systems at small frequencies. 

We begin by making the ansatze 

'ip±{u) = f{u)^ J^±{u), 

\ m J (3.1) 

jr^(M) = F^\u) + uF2\u) + u^F^\u) + 0{u^), 

for the gauge invariant fields, where we are using the dimensionless variables 

To ^ u ^ m ^ fi 

u = —, oj = —, m = —, fi = —. (3.2) 

r ro ro Tq 

We have factored out an oscillating function that corresponds to imposing ingoing 
boundary conditions at the black brane horizon, as well as an overall w-dependent 
function such that the leading terms F^\u) will be independent of u. To determine the 
functions F^\u)^ we substitute the ansatze (3.1) into the equations of motion (2.4) and 
expand as a power series in a). We then solve order-by-order for the functions 


- 12 - 




demanding that (m) is regular and equal to a frequency-independent constant at the 
black brane horizon u = 1. At leading order in a), we hnd that F^^\u) = Ci is a 
constant which we will set to 1 for convenience. At higher orders, are non-trivial 
functions of u that satisfy equations of the form 

Ff{u) + g^^\u)Fi^\u) = S^^\u). (3.3) 




u), with the source 


That is, they are hrst order, linear inhomogeneous equations for F± 
terms (u) depending on the solutions at lower orders in the perturbative expansion. 
In principle, exact integral solutions to these equations can be found by using the 
method of integrating factors. We present the technical details of the perturbative 
solutions in appendix B, and focus on the physical consequences in the following. 

Using (2.14), the diagonal correlators are 


{J±J±){u:) = 


(74 


7-) 1 + uF^\0) + u^F^\0) + 0(cj3) ^ 


(3.4) 


and these can easily be combined to give the full conductivities using (2.13). The DC 

f 1)^ 

conductivities are controlled only by F^^> 

(0), which can be analytically determined 
exactly as a function of (see appendix B). Using these expressions, we recover 

the results (1.6) for the DC conductivities which were discussed extensively in section 
2.3. This provides an independent check of the horizon formula results of [14, 18] from 
a Kubo formula calculation. 

To understand the physical origin of each contribution to the conductivities, we 
must determine their frequency dependence by working to higher orders in the pertur¬ 
bative expansion u. We were not able to obtain analytic results for general T, p, m at 
these higher orders.^ Instead, we have focused on the limit of slow momentum relax¬ 
ation m 1, and computed the conductivities in a perturbative expansion at small 
u and rfi^, with a) ~ -C 1, i.e. at long timescales, comparable to the momentum 
lifetime. For some of the coefficients of the terms in (3.4), we have only been able to 
obtain analytic answers perturbatively in p. However, our hnal result for the coherent 
and incoherent contributions to the DC conductivities will be non-perturbative in p. 
The details of the perturbative solutions are given in appendix B. Note that we are 
working with dimensionless variables normalised by tq: for m 1, ro ~ T when /i ^ T 
and ro ~ /i when T ^ /i. 


3.1 The incoherent contribution 

In this subsection, we will focus on the (J+J+)(a;) correlator at small m. Recall that 
in the strict m = 0 limit, J+ is the current (2.15) which decouples from momentum 

note that exact results in m, T can be obtained for the neutral ^ = 0 state. 
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and which is therefore completely incoherent. Our perturbative calculation, described 
in appendix B, yields the J+ conductivity® 


+ - 




u 


[ctq + 0(m^)] + [/3i + 0(m2)] u + 0{u^) 

1 + [/S 2 + 0{rh‘^)] a; + [/S 3 + 0{rh‘^)] a;^ + 0{oj^) ’ 


(3.5) 


where 

. ^2 •~4 

A = (v^tt + giogsj “ 216 “ 63 log 3^ + ^ (jVSti - 72 log3j + 0(/i®), 

/32 = (^^/Svr + 9 log 3j ^ (^igVSTr - 9 log3j - - 41og3j +0{fl^), 

/^s = -^ (tt^ + 6y37rlog3 - 27 (log3)^j + 0(/i^). 

(3.6) 

The nature of the transport of the current J+ is encoded in its pole structure: coherent 
transport is caused by a parametrically long-lived excitation. In our system, this would 
be a Drude-like excitation due to the slow relaxation of momentum, which has the 
dispersion relation O ~ [7, 21]. 

It is clear from (3.5) that there is no such long lived excitation transporting the 
current J+. Our calculation shows that the longest lived collective excitations in this 
sector have microscopic lifetimes ~ Therefore the contributions of (J+J+) to the 
thermoelectric conductivities (2.13) are all incoherent. To the order of the perturbative 
expansion to which we are working, the contribution of these incoherent processes to 
the full thermoelectric conductivities of the system is 


cr+(w) =aQ + 0(a), m^), a+(a;) = -^ctq 0(a), m^), 

-+/ \ ~ 2\ (12 —/i^)^ 

^■*■( 0 ;) = —(To -|- 0(a;, m ^), an = - 

K ) T ^ l Q 9(4 + /12)2 


(3.7) 


This is entirely in agreement with the incoherent contributions to the conductivities 
predicted in (1.3) by the hydrodynamic effective theory of [1]. 


3.2 The coherent contribution 

We now turn to the (J_J_)(a;) correlator. We know that this must have a coher¬ 
ent component at small m since the conductivities themselves do. However it is not 

®We normalise by a factor of due to the ubiquitous appearance of such a factor in (2.13) at 
small m. 

®We cannot give quantitative results for the lifetime, as frequencies a) ~ 1 are outside of the range 
of validity of our perturbative calculation. 
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clear, a priori, if J_ is transported completely coherently (to the order in to which 
we are working), or whether it has both coherent and incoherent components. Onr 
pertnrbative calcnlation, described in appendix B, yields the J_ condnctivity' 


oj 


where 


[ai + a2rfi^ + 0{m^)] — [bi + 0(m^)] iu + 0(0'^) 

fn^ — [ci + C 2 fn^ + 0{fri‘^)] iO + [di + 0{m'^)] a)^ + 0(0^) ’ 

(3.8) 


CLi — 

8/x^ (12 + /x^) 

9(4 + p2)2 > 

^ (\/37r - 18 + Qlogsj + ^ (27 - 4^7r - blogsj + 0(/i®), 

= 4 (4 + ) 

^ ~2 ~4 

C 2 = (9 + VStt “ ^logsj + ^ (72 - ^ (a/Stt - 84 + 3log3j + 0(/i®), 

* ” “ 6+ “ '^) + M “ ‘“®^) + 

(3.9) 

Note that we are working to the hrst snbleading order in the small Cj ~ m? expansion 
in both the nnmerator and the denominator. 

As expected, this correlator has a pole at cD ~ —ifh?. This pole corresponds to the 
existence of a parametrically long-lived Drude-like excitation dne to the slow relaxation 
of momentnm, and will give coherent contribntions to the thermoelectric condnctivi- 
ties. Onr pertnrbative calcnlation allows ns to determine snbleading corrections to the 
location of this Drnde-like pole 00 d 


ud = 


3 (4 + /i2 


-m 


— (9 + VStt - 9 log3) + ^ (-198 - x/Svr + 81 log 3 ) 




2592 


187 — Ga/Stt — 54 log3j + 0(/i®) + 0(m®) 


(3.10) 

The O(m^) term agrees with that of [21], as it shonld (given the similarities between 
the gravity theories nnder stndy here and there [14]). The ft = 0 limit of this term is 
also in agreement with [7]. At 0(?Ti^), we could only calculate the location of the pole 
perturbatively in ft, as is clear from the result (3.10). As a check of our calculation, a 


^See footnote 5. 
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Figure 1. A comparison of our analytic result (3.10) for the O(m^) correction to the location 


of the Drude-like pole (blue line), and the exact location obtained numerically for m = 1/10 
(black dots). There is excellent agreement at small jl where our perturbative analytic result 
should be accurate. 

comparison of this analytic result with a numerical calculation of the pole location is 
shown in figure 1: there is excellent agreement for appropriately small values of fi. For 
the neutral /i = 0 case, interestingly, the location of the pole to is equivalent to 

that of the transverse momentum diffusion pole in the translationally invariant system 
dual to the Schwarzschild-AdS 4 black brane (see table III of [31]), after replacing m 
with the wavenumber q. 

From this pole, we can identify the momentum relaxation rate F as F = iroCon, 
which can be written (perhaps more illuminatingly) as 



y/Sn — 9 log 3 9/i^ (log 3 — 2) 

%7iV^ 256?r*r-> 



(3.11) 


where the thermodynamic quantities are m-dependent. There will, of course, be other 
poles in the correlator with decay rates ~ but our perturbative calculation is not 
able to accurately capture these. 

The current J_ clearly has a coherent contribution to its transport, due to the 
existence of the Drude-like pole. Generically, we would also expect there to be an 
incoherent component to its transport. To quantify this, we calculate the residue of 
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the Drude-like pole in the conductivity of J_: 


= lim {u — u;d)T, (o;). 


(3.12) 


The natural dehnition of the contribution of the coherent excitation to is then 
— ZdICjd- For a conductivity of the form (3.8), this yields 



(3.13) 


as the coherent contribution to Sbc- At order all of the DC conductivity comes 
from the Drude peak, but at the hrst subleading order, 0(m°), this is not necessarily 
the case. The term in brackets in (3.13) indicates a part of which does not come 
from the Drude peak i.e. it is an incoherent contribution. However, substituting in the 
explicit expressions for ai,6i,ci and di for our system (3.9), this potential incoherent 
component of the DC conductivity vanishes identically! This indicates that, to the 
order in p to which our result (3.9) is valid, the entire DC conductivity of J_ comes 
from the Drude-like excitation, up to and including the hrst subleading order in the 
fh? expansion. 

In fact, we can show that this is true to all orders in fi. Although we do not know 
individually how ai,6i,Ci and di depend upon fl, it is easy to check that the precise 
combination appearing in the brackets in equation (3.13) must vanish, by demanding 
that in the strict m = 0 limit, we reproduce the hydrodynamic results of [26]. This 
assumption of continuity of the hydrodynamic limit is manifestly true up to 0(/i^), 
and we believe it should be true to all orders. In summary: up to and including the 
hrst subleading order in the expansion, the entire DC conductivity of J_ comes 
from the Drude-like excitation - there is no incoherent component at this order. At 
higher orders in the expansion, we expect to be a sum of both coherent and 
incoherent contributions. 

The contributions of this sector to each conductivity are therefore totally coherent, 
to this order, and given by 



(3.14) 








where ctq is given in (1.7) and F is given in (3.11). At higher orders in the small 
a), expansion, we expect that the contributions of J_ to the conductivities will be 
comprised of both coherent and incoherent pieces. 


3.3 Discussion 


Collecting the results (3.7) and (3.14), the thermoelectric conductivities are given by 
equations ( 1 . 8 ) in the limit of slow momentum relaxation. As we have demonstrated, 
the coherent part of each conductivity comes solely from J_ at this order, while the 
incoherent part comes only from J_|_. For an easier comparison with the memory matrix 
and hydrodynamic formulae (1.2) and (1.3), we can change variables from to F and 
write the conductivities to subleading order in a small a; ~ F expansion 


a (a;) 
a (cn) 


K (cn) 


^ + F(l-aQ + A/x^) + 0(F^n;F,n;^) 

Y — ibj 

^ + r (gJg + 47Tn\) + 0 {T\ojT,w^) 

Y — iuj 


+ <Jq + 0(u, F), 

- ^(^Q + 0(u;,Y), 


+ F + dTTsTA) + 0(F2, ujY, oj^) 


(3.15) 


In the translationally invariant limit F = 0, these agree with the results of [26], and the 
hydrodynamic formulae (1.3). When F 7 ^ 0, they agree with the memory matrix results 
(1.2) and hydrodynamic results (1.3) at leading order in the small cn, F expansion, 
but not at subleading order. The subleading corrections in (3.15) are comprised of 
two independent pieces: an incoherent contribution, and a coherent contribution (a 
correction to the weight of the Drude peak). The hydrodynamic results (1.3) correctly 
capture the incoherent contribution but not the correction to the Drude peak. The 
memory matrix results ( 1 . 2 ) can also be extended to include an incoherent contribution 
[ 8 ], but not yet the correction to the Drude peak. Since the subleading correction to 
the Drude peak enters at the same order (in a small F or small m expansion) as the 
incoherent contribution in our holographic theory, it is important that these effective 
theories are extended to incorporate this correction to the Drude peak. 

In the limit of zero chemical potential, the conductivities are given by 


a{uj) 
a (cn) 

k (cn) 


1 , 

0 , 

^ + dvrsTAF 

e+p_ 

Y — iu 


(3.16) 
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to this order. The electric conductivity a is totally incoherent since J decouples from 
P when /i = 0, while the thermal conductivity R is totally coherent at this order, 
conhrming further the results of [7]. 

In the introduction, we noted that previous works have tried to identify the m- 
independent contribution to anc (which numerically is equal to 1) as being the inco¬ 
herent component of the electrical conductivity. As is clear from our results, this is not 
the case. However, the m-independent contribution in this theory can be identihed as 
being the DC value of the electrical current in the absence of heat flow [18] 

o'dc\q=o = ctdc -—^ = 1. (3-17) 

kdc 

This result can be generalised to more complicated holographic theories in a natural 
way [18, 32]. With our results, we can revisit this computation and determine how this 
conductivity depends on frequency, hnding 

2 

-aQ + 0{u,T), 

(3.18) 

= 1 - 1-0 - -—,.. .^ , 

\ T — lu J 


^M|q=0 = 


lUJ 


(e + p)' 

g2T2 


(Tq + 0{u:,T) 


+ 


(e + p) 

g2T2 


to the order to which our calculations are valid and recalling the value of ctq (1.7). This 
conductivity is totally incoherent to this order, and does not have any contributions 
from subleading corrections to the Drude peak. It would be interesting to determine 
whether this is also the case to higher order in the expansion, or whether 
has contributions ~ T^/(T — iu) etc. Similarly, the heat conductivity in the absence of 
electrical current is 

2 

-ctq -|- O {u, T) , 

(3.19) 




T 


T — iu 
s^T 


Ts^ {e + pY 






-|- O ( T, w. 


n^T 

r2 

V-iu 


+ O (w, T) 


+ 


(g + p) 

n^T 


which is totally incoherent to the order to which we are working. We note that the 
absence of any leading order contribution ~ r°/ (T — iuj) to these conductivities is as 
expected from [11]. 


4 Outlook 

We have shown that the transport of heat and charge in the state with momentum 
relaxation, dual to (1.5), can naturally be expressed in terms of the two currents J±, 
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given in (2.12), which diagonalise the thermoelectric conductivity matrix (1.1) for all 
values of m. In certain limits, the form of these currents can be used to understand 
the physical processes underpinning the transport properties. In the limit of very fast 
momentum relaxation [m —)■ cxd), the heat and electrical currents decouple, as they do 
in the charge conjugation symmetric limit. In the limit of no momentum relaxation 
(m = 0), the decoupled currents are the coherent energy current Je (which is equal to 
the total momentum P), and the current Je — ^J which is completely incoherent, as 
it decouples from the total momentum P [27]. 

We have analytically computed the low frequency behaviour of the conductivities 
in the limit of slow momentum relaxation (small m). In this limit, the decoupled 
currents J± are still controlled by qualitatively different physical processes. To the first 
subleading order at small m, J_ remains coherent, i.e. it is controlled by the momentum 
relaxation timescale of the system, while J+ remains incoherent, i.e. it is controlled by 
the intrinsic relaxation timescale of the system. There is a smooth m —?■ 0 limit. 
The two independent contributions combine in a very non-trivial way to form the DC 
conductivities (1.6) — it is not easy to guess how the DC formulae should be divided 
up into coherent and incoherent contributions without any other information. 

Our results highlight the fact that subleading corrections to the Drude weight enter 
at the same order (in m) as the leading incoherent contribution to each thermoelectric 
conductivity. The apparent discrepancies between the holographic DC conductivities 
and those of the memory matrix or hydrodynamic descriptions are due to the neglection 
of corrections to the Drude weight in these effective theories. 

There are several directions which are worth pursuing further: 

Spatially resolved transport We have considered the transport of the spatially 
uniform components of the charges and currents. A natural extension would be to study 
the transport of the non-zero wavenumber k harmonics, to understand how charge is 
transported over different distance scales. In the limit of slow momentum relaxation, 
we expect that, at low energies, J_ will be transported by sound at short distances 
(large k) and diffusion at long distances (small k), as was observed in [7] for a zero 
density system. In contrast to this, we expect that J+ will be transported by diffusion 
at all distance scales, due to its incoherent nature. 

Magnetotransport Building on [1], a number of recent articles have revisited the 
problem of magnetotransport with momentum relaxation by computing the thermoelec¬ 
tric and Hall conductivities either holographically [5, 33-35] or with memory matrices 
[8]. To resolve the discrepancy between the hydrodynamic, memory matrix, and holo¬ 
graphic DC calculations, it would be worthwhile to adapt our techniques to calculate 
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the frequency dependent conductivities at non-zero B. Extending our calculations to 
non-zero B would also allow us to examine whether the Hall angle receives contributions 
from both coherent and incoherent processes and how this relates to the interpretation 
of its temperature scaling in terms of two timescales [5]. 

It was recently proposed [36] that the timescale setting the resistivity scaling of 
the strange metallic region of a certain iron pnictide compound is proportional to the 
square root of a sum of squares of the temperature and magnetic held. With this in 
mind, it would be very interesting to determine the dependence of the appropriate 
timescale (momentum relaxation rate or diffusion constant) on the magnetic held in 
holographic systems. This could be done by adapting the methods we have used here. 

More general theories It would be very worthwhile to extend our work to more 
general holographic theories with slow momentum relaxation, in which a hydrodynamic 
limit exists at non-zero temperatures. This should be the case when a Drude-like pole 
dominates the correlators at sufficiently low energy scales. Holographic theories can 
exhibit branch cut formation in the T —)■ 0 limit, due to a coalescense of poles with 
decay rates differing by ~ T. Although our analysis will not capture these poles, a 
hydrodynamic limit should be valid when cj, T <C T, as this is when the Drude-like 
excitation is parametrically longer lived than the rest. 

In theories where there is a neutral scalar which can run logarithmically in the 
interior of the geometry [16, 17, 19], we would expect that the temperature scalings of 
the coherent and incoherent contributions to the conductivities can be different from 
one another. With this additional hierarchy of scales, it may be possible to hnd states 
with slow momentum relaxation where the effects of the incoherent contribution are 
parametrically larger (or smaller) than corrections to the Drude peak. The method 
used in [22] may be useful for more general theories. 

Another question is the sensitivity of our results to the choice of momentum relax¬ 
ation mechanism: would they be modihed if we had instead used random-held disorder 
[3, 4, 22, 37-39], or homogeneous [19, 40, 41] or inhomogeneous lattices [20, 32, 42] to 
break translational invariance? Furthermore, if we had broken translational invariance 
with electrically charged, rather than neutral, operators, would this affect the nature 
of transport in the system? In particular, we interpreted the decoupling of J and Q 
at large m as being a consequence of the state’s thermodynamics becoming dominated 
by the neutral scalar degrees of freedom. Does the same decoupling occur (at low 
frequencies) when these neutral operators are not present? 

A qualitatively different class of holographic systems with hnite conductivities are 
probe brane systems, whose DC electrical conductivity can be written as the square 
root of the sum of two terms [43], one of which is often interpreted as a ‘Drude-like’ term 
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(and can also be compnted from the drag force on the charge carriers), the other as a 
‘pair creation’ term. It wonld be interesting to verify this interpretation by analytically 
compnting the low freqnency, linear response condnctivity in snch a system, as we have 
done here. 

Effective theories of thermoelectric transport with slow momentum relax¬ 
ation Onr computation has highlighted what needs to be done to rehne existing 
effective hydrodynamic [1] or memory matrix [8] theories of transport in the presence 
of slow momentum relaxation, such that they are consistent with the holographic com¬ 
putations of DC conductivities. These effective theories should be extended to take into 
account order T corrections to the weight of the Drude peak. These produce 0 (r°) 
corrections to the DC conductivities, which are the same order as the incoherent ctq 
contributions. This is an excellent example of how gauge/gravity duality can contribute 
to the understanding of transport in strongly correlated systems in general, by provid¬ 
ing a consistent and reliable framework from which effective theories can be extracted, 
or to which effective theories can be compared. 
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A The on-shell action 


The on-shell action of the theory (1.4), to qnadratic order in the perturbations (2.2) 
around the solution (1.5), is 


s = J mhf^\u) + iux2^^\uj) + ^ay^^\-u)ay^^\u) 

(A.l) 

where we have expanded a generic held perturbation 6ip{r,uj) near the boundary as 


6(p{r,u) = 


(A.2) 


and set the scalar operator source term x^\^) fo zero. From this, we can use the 
standard AdS/CFT dictionary [44] to calculate expressions for the retarded Green’s 
functions of the operators dual to each held perturbation, in terms of a® (cn), a|,^^(a;), 
etc. These can then be rewritten in terms of the near-boundary expansions of the 
decoupled variables using their dehnitions (2.5). To compute the subtracted correlators 
(2.10) that enter in the Kubo formulae (2.1) for the conductivities, we must subtract the 
retarded Green’s functions when a; = 0 and fc —)■ 0, where k is the wavenumber of the 
perturbation in the y-direction. These were obtained by computing the on-shell action 
for huctuations of this kind, yielding the expressions (2.9) for the subtracted correlators. 
A non-trivial consistency check of our calculations (including contact terms) is that, 
after solving the equations of motion and substituting these solutions into (2.9), we hud 
that the conductivities are free of i/u poles, as should be the case on physical grounds. 


B Details of the perturbative calculations 


In this appendix, we give details of the perturbative solutions for the functions fF±{u), 
dehned in (3.1). At leading order in a), the solutions which obey the correct boundary 
conditions at the horizon are simply constants F^\u) = C±. The value of these 
constants is unimportant and will cancel out in the hnal answers for the conductivities, 
and so for convenience we set C± = 1. 

At 0(a)) in the expansion, we can formally write the solutions as integrals 

4 (m 7 ± -I- p)^ (p^ -I- 2m^ — 12) -|- 4xh±(x) [nF (6a; — 4) -f a; {p? (4a; — 3) — 12)] 


F^^\u) = / dx 


i{x — 1) h±{x) {jF + 2m? — 12) [—4 — 4x -|- a;^ (—4 -1- 2m‘^ + xjj?)] 

(B.l) 
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where 


h±{x) = (m 7 ± + . (B.2) 

j" 1 \ / 

From these integrals, it is straightforward to analytically calculate the constants F± (0) 
that control the DC conductivities 



.(m7± +/i)^ 
^ ~ 2 2 


(B.3) 


However, we could not do the integrals analytically and hnd exact expressions for 
7^7)- We are primarily interested in the small m limit of the conductivities, and 
thus the small fh limit of the integrals. For 7^7)) if is straightforward to expand the 
integrand at small m, and hnd that the leading order term is of order This means 
that there is no Drude-like excitation in the conductivity of J+, which is therefore 
incoherent. For our purposes, this is all we need to know. For completeness, we note 
that it is not possible to integrate the leading term of the integrand analytically for 
general p, but that it is possible in a small ft expansion; 


F<‘'(0) = 


• ~ 2 

^ (VSti + 9 log s) + ^ (^19A/37r - 9 log3j + 0(7) 


+ 0{m^ 


(B.4) 

For 7^7)) if is more complicated. The small m limit of the integrand is singular 
due to the form of the function h_{x) ~ (m^ + in the denominator: the limits 
rfi —)■ 0 and a: —)■ 0 do not commute. If we hrst send —)■ 0, the integrand diverges 
when X —)■ 0. To correctly evaluate the small m limit, we must take it small but non¬ 
zero, so that we accurately include the contribution from integrating over the region 
0 < X < To do this, we change the integration variable to x = rh'^y before 

expanding the integrand at small rh and integrating the leading term in this expansion 
to give 

3i (4 -|- 7'' 


7^7) = 


4m^ 


-I- O(m^). 


(B.6) 


This change of variables is only useful for giving us the leading term: it does not allow 
us to accurately extract any of the subleading terms in m. We have checked that this 
technique is reliable by explicitly doing the integral numerically and comparing it to 
our result (B.5). The consistency between our analytic pole location and the exact one 
determined numerically (see hgure 1) is also a check of this. The hrst correction to 
(B.5) is given below in (B.6). 

At second order in the small O expansion, things are even more complicated and we 
can only get analytic results for (m) by performing a double expansion at small m 
and small jl. The strategy is as follows: we expand the integrand of (B.l) to the second 


24 






subleading order in /r, and integrate each coefficient to obtain an expression for F±\u) 
which is perturbative in ft and exact in fh. This enters as a source in the equations 

/o\ / 

of motion for Ff (u), for which we can write down formal integral solutions which 
are much too lengthy to include here. We then expand these integrands to the same 

~ f 2^) ^ 

order in p and again integrate term-by-term to obtain expressions for {u) which 
are exact in m but perturbative in ft. The hnal step is to integrate these expressions, 
but we could not do this analytically, even in the small ft expansion. Since we are only 
interested in the leading order behaviour at small m, we expanded each term in the 
small /i expansion of the integrand to the lowest order in m. For F_ (u), this again 
was preceded by a rescaling of the integration variable x = fh?y due to the singularity 
of the m —)■ 0 limit of the integrand. The results are as follows 


(0) = — (ySn + 9 log 3j ^ (bVSTt - 63 log 3j -F 0{fi^ 






F^\0) = 


^ft^ 


^ (VStt- 18 + 91og3) - ^ (9 + 2\/37r- 121og3) + 0{ft 


24 


-h 0{m ), 


216 


TT^ -|- gVStt log 3-1-27 (log 3)' 


(TT^ — 9 (log 3)^ -I- 18-\/37rlog3 -1- ^ ( ^ ) ) + 


432 
F?1(0) = 


^ (l8-y37r-91og3) (6 - 2y37r + 2log3) + 0(/i' 
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-F O(m^), 


-I- 0(m° 


where is the polygamma function. A byproduct of this analysis is that we 

obtain the 0(m°) correction to (B.5), perturbatively in ft: 


f[^\o) = 


3i (4 -|- /i^ 

4^2 


+ 


— 9 + 
18 


y37r-91og3) + ^ (^y37r-72-91og3 




- f-84 a/Stt 3 log 3) -h 






(B.6) 

The calculation we have described is quite complex and involves taking two limits 
(small ft and small m) which, in principle, may not commute. But there are a number 
of consistency checks we have performed to make sure the expressions above are correct. 
For example, the location of the Drude-like pole (3.10) is sensitive to the value of f 1^1(0), 
the hnal quantity derived in the procedure above, and our analytic expression agrees 
with the exact numerical result (see hgure 1). The form of the conductivities in the 
m = 0 limit also depend non-trivially on these coefficients (as described in section 3.2), 
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and we have checked that we recover the correct results in this limit. Finally, where 

possible we have numerically computed the integrals and checked that the results are 

consistent with our analytic expressions. 
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